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Abstract 

We investigate the stability of circular material orbits in the analytic galactic 
metric recently derived by Harko et al. (2014). It turns out that stability 
depends more strongly on the dark matter central density po than on other 
parameters of the solution. This property then yields an upper limit on po 
for each individual galaxy, which we call here such that stable circular 

orbits are possible only when the constraint po < is satished. This is 

our new result. To approximately quantify the upper limit, we consider as a 
familiar example our Milky Way galaxy that has a projected dark matter radius 
i?DM ~ 180 kpc and hnd that pg^^®"^ ~ 2.37 x 10^^ M 0 kpc“^. This limit turns 
out to be about four orders of magnitude larger than the latest data on central 
density pg arising from the fit to the Navarro-Frenk-White (NFW) and Burkert 
density profiles. Such consistency indicates that the EiBI solution could qualify 
as yet another viable alternative model for dark matter. 
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The issue of dark matter is at the centerstage of modern astrophysics. Sev¬ 
eral theoretical models and cosmological scenarios for dark matter exist in the 
literature[l], among which the possibility of perfect fluid dark matter within the 
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framework of general relativity has already been explored in the literature [2,3]. 
A similar possibility has been recently investigated within the framework of the 
Eddington inspired Born-Infeld (EiBI) theory by Harko et al. [4]. 

Using a tangential velocity profile [5] giving flat rotation curves at large 
distances and setting the cosmological constant A to zero, they obtained, in the 
Newtonian approximation, a new galactic metric and theoretically explored its 
gravitational properties including the dark matter density distribution. 

However, the numerical values of the crucial parameter k (denoted hy k = 
2 -Rdm/'^^) or equivalently the dark matter radius i?DM cannot be determined 
from the theory alone — it has to be supplied by the observed data. This i?DM 
is expected to make its appearance also in the stability of orbits. 

The purpose of the present Letter is to investigate the stability of circular 
material orbits in the EiBI galactic metric [4], which is by no means obvious from 
the metric itself. It turns out that the stability is exclusively sensitive to the 
variation of the dark matter central density po leading to an upper limit, which 
we call here such that stable circular orbits in the EiBI are possible only 

when the constraint po < is satisfied. We shall illustrate the inequality 

by the example of our own Milky Way galaxy, which shows that pg^^®"^ is about 
four orders of magnitude larger than the latest data on pg arising from the fit 
to NEW or Burkert profiles. We shall take units such that G = I, c = I, unless 
otherwise specified. 

The salient features of the EiBI dark matter model are as follows: The action 
is 

•SEiBI = ^— [ d'^x 
Sttk J 


\J- IPpi/ + SttkA^j.] - Av^ 




( 1 ) 


where A is a dimensionless parameter and «: is a parameter with inverse di¬ 
mension to that of the cosmological constant A. In the limit «; —)■ 0, the 
Hilbert-Einstein action is recovered with A = SttkA -|- I, where A is the cos¬ 
mological constant. Harko et al. [4] developed spherically symmetric solution 
assuming A=l^A = 0so that n can have nonzero values. The descrip¬ 
tion of the physical behavior of various cosmological and stellar scenarios was 
assumed to be controlled by the single parameter k. The galactic halo is as¬ 
sumed to be filled with perfect fluid dark matter with energy-momentum tensor 
Tr-’' = pg^'^ + (p + p)U^U''^ g^iuU^U'' = — 1. There are two metrics in the EiBI 
theory, the physical metric g^^^ and the auxiliary metric and the tangential 
velocity profile is taken as [5] 
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where r is the standard radial coordinate (as defined by Eq.(9) in [4]), ropt is 
the optical radius containing 83% of the galactic luminosity. The parameter rp 
is defined as the ratio of the halo core radius and ropt) arid Uoo is the asymptotic 
constant velocity. Under the Newtonian approximations that the pressure p ~ 0, 
Stt^p <C 1, and (r/ropt)^ ^ rg, the EiBI field equations yield the Lane-Emden 
equation with polytropic index n = I, which has an exact nonsingular solution 
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for dark matter density as [4] 


p(^\r) = K 



( 3 ) 


where (0) = K = po is the central densityj^ Assuming that the halo has a 
sharp boundary -Rdm, where the density vanishes such that p^^^Rbm) = 0, one 
has 

( 4 ) 


( 5 ) 


^DM = TJ" A / 2 • 

Thus, the mass profile of the dark matter is 

4^3, 


/•r Af>3 

M{r) = An p^^\r)r‘^dr = -[sin(r) — f cos(f)] , 

Jo 


where the dimensionless quantity r = 7rr/i?DM- The average velocity dispersion 
of dark matter particles in the constant velocity region in the present model is 

[ 4 ] 


0-2 = = 


'tg 


1 - ^73cCi(r) 

sin(rj 


( 6 ) 


where Ci(z) = — /. 


_ COO cos(t) 


dt is the cosine integral function. 


The approximate physical metric has been derived as [4] 

dr^ = -B{f)dt^ + A{f)df^ +f^C{r){d0^ +sm'^ 0d(j)^), (7) 


A{r) 


B{r) = 

= 6^0 

f Rdm \ -2 , 2 
r +rg 
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1 — ^ sin(r) + pg cos(r) 

1 — ^ sin(f) 


C{f) 



1 - ^ sin(r) 
r 


( 8 ) 

( 9 ) 

( 10 ) 


where is an arbitrary constant of integration (which we set to unity) and 
the dimensionless quantity pg = ^ jsjote that the surface area of a sphere 

at the boundary of dark matter halo dehned by f = tt, has the value S = 
Anr^Cif) = Anr^ = 47ri?p]y[, which is just the spherical surface area in 

’’standard coordinates”. Thus the dark matter radius i?DM can be identified 
with standard coordinate radius. 

^The solution (3) of the Lane-Emden equation and its connection to dark matter were 
first pointed out in Refs.[26-28]. The same density profile is used also in the context of the 
Bose-Einstein Condensate (BEC) simulation of dark matter [29]. 
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To analyze the stability of circular orbits, one needs to analyze the second 
order derivative of the concerned potential, which we wish to do here. To find the 
potential V, note that the four velocity C/“ = of a test particle of rest mass 
Wq moving in the halo (restricting ourselves to 9 = 7r/2) follows the equation 
guaU'^U'^ = —TOq that can be cast into a Newtonian form in the dimensionless 
radial variable r (= 7rr/i?DM) as 


dr\‘ 


Tr) 


which gives, for the metric Eqs.(7)-(10), the potential 
V{r) = 


( 11 ) 

( 12 ) 

(13) 


where the constants E and L, respectively, are the conserved relativistic energy 
and angular momentum per unit mass of the test particle. Circular orbits at 
any arbitrary radius are defined by r = i? = constant, so that ^ \r='R= 0 
and, additionally, ^ lr=S= 0- From these two conditions follow the conserved 
quantities: 

= ^ (14) 
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and using it in V (R) = —E"^, we get 


where 


^ -z’ 


X = —k^R v^{R — Pq sin R)^ 


(15) 


(16) 


Y=(kR^ + 2rgr2pt 


—2 \ ’'o 

kR \ 


° 2r2 

^'opt 


(/9gi?COSi? + Po sini? — 2i?) (17) 


Z = ^kR^ {l-2vl^)Po sin R+(^kR^+ PqR cos R 

-ARrlrlp^ - 2 kR^ (l - . (18) 

Putting the expressions for and E'^ in Eq.(12), we find the complete 
expression for V. The orbits will be stable if V” = \r=R< 0 ^nd unstable 
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if V" > Q. The expression for V" is 




2v'^ [R + PqR cos R — Pq sin R) 

R^ (^kR^ + 2r2r2pt^ Z 

32i?V^r2pj + SkR'^ (1 - + hp^rgr^p^ (^1 + R^'j 


+kpIR^ (l - 2^;^ + 3i?") - 2i?" {wr^r^p^ + nR^ (5 - 2r;^ 
+ {2 - 3) r^r-^pt - kR^ (1 - 2 r;^ - R^'j | pi cos ( 2 i?) 

- {36r^r^pt + 4:R^rlrl^^ + 6 kR^ + 2kR^ - 12Ki?\^| pgi? 

+ (Sj’o^opt + Po^sin (27?) . 



sini? 


(19) 


From the above expression, it is absurd to straightforwardly draw any con¬ 
clusion about stability or otherwise of the circular orbits. Clearly, much will 
depend on the parameter ranges chosen on the basis of physical considerations. 
While other parameters can be reasonably assigned, the as yet unknown param¬ 
eters are the dark matter radius k (= 2R^y^/ tt^) and the dimensionless central 
density pg (= 8poi?Qi^/7r), again depending only on k. In the first order ap¬ 
proximation, the density distribution in the dark matter has been assumed in 
[4] to be low such that ?>ttGkp ''^^/<< 1, but the central density pg could still 
be large since |sin(a;)/x| < 1 [see Eq.(3)]. The question therefore is how large 
or small could it be, or turning around, could there be any upper limit on pg 
imposed by the stability criterion? 

The answer to this question is yes and can be found graphically. We find 
that V" is indeed very sensitive to changes in pg leading to different upper limits 
PgPPer for different galactic samples such that stable circular orbits are possible 
only when pg < pg'’^®’^ in that sample. Different pg^*’®’^ results from the fact 
that i?DM changes from sample to sample, as it should, and thereby leads to 
different (though not too different) values for k. For illustrative purposes, let us 
fix typical values for a galactic sample, say UGC 0128, a low surface brightness 
galaxy of moderate size, with the last observed scattering radius occurring at 
.Riast = 54.8 kpc. Since direct observational data on i?DM is yet unavailable 
for any sample, we approximately determine it by projecting beyond i?iast the 
observed tendency of continuous decline of the velocity dispersion to the zero 
value (see Fig. 6 of Ref. [7]). For the present sample, one can then read off a value 
i?DM ~ 88 kpc, which corresponds to k = 1.57 x 10^ kpc^ = 1.49 x 10'^® cm^ 
(conversion: 1 kpc= 3.085 x 10^^ cm). The other relevant parameters within the 
Newtonian approximation ar^ rg = 0.044, Uoo = 0.000001, ropt = 27.6 kpc, 

^The range of r and ropt is chosen so as to ensure r/ropt rg. The usual formula for rg 
evaluates to 0.044 for the sample UGC 0128, while ropt — 4i?o- For this sample, Rq = 6.9 kpc 
[7] such that ropt = 27.6 kpc. We emphasize that exact values of these parameters including 
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the dimensionless radius R (= ttR/Rdm) is chosen in the range R G [ 0 . 57 r, 7 r] 
and the dimensionless density parameter in the range Pq G [0.257r, O.Stt]. 

Graphical analysis reveals a remarkable result: While V remains practi¬ 
cally insensitive to the variation of the parameters («:,ro, v^o, ’"opt) within the 
periphery of Newtonian approximation, it is greatly sensitive to the variation 
of the remaining parameter po- Figs.l and 2 respectively show that, for val¬ 
ues of Pq > 0.94, there is instability in the entire or partial range of the halo 
radii R, while Fig.3 tells us that there is an upper limit occurring at Po^^***^ = 
0.94 = A^pp^’^tt, where A'^pp®"' = 0.299, su^h that for p^ < p^PP®’' all circular 
orbits in the entire chosen radial range for R are stable. Note that this value of 
yupper game under the change of parameter values (k, tq, Voo, ’"opt) 

in a given single sample, and so PqPP®^ is quite robust for that sample. 

Rewriting in terms of po, we have 


upper _ 

Po — 


P^PP^TT 

8Rdm 


yupper 

4k 


( 20 ) 


which implies that kpOpp*^*^ = constant, that is, the larger the dark matter radius 
Rtom, the lower the value of PqPP®^. This is an interesting feature of the EiBI 
model. Plugging the values of A^pp®*^ and k, we find that the constraint pg < 
—upper jjg^mediately translates into an upper limit on pg such that for 


Po < Po 


upper 


( 21 ) 


all circular orbits in the chosen range for R are stable. We have verified that 
-upper jjggpg |.p ^j^g ganig value (pg^^™ = 0.94) for many samples listed in [6,7], 
the only difference to pgPP®"^ thus comes from the various Rdm characteristic 
of various samples, as expressed in Eq.(20). Eor UGC 0128, pgPP®*^ ~ 9.9 x 
10^^ M 0 kpc“^, noting the conversion: 1 cm“^ = 1.98 x 10^® M 0 kpc“^. Thus 
the prediction is that, as long as po of any galaxy obeys the stability induced 
constraint ( 21 ), the circular material orbits in the halo around such galaxies 
will be stable up to a maximum radius R = Rum- Otherwise, the orbits will 
be unstable so that there will be no dark matter around galaxies that have po 
exceeding the ppPP*^*^. Eor further support to the constraint, we may consider a 
large LSB sample U11748 with a projected dark matter radius Rum ~ 106.18 
kpc, which leads to ppPP®*^ = 6.81 x 10^^ M 0 kpc“®. This upper limit is quite 
comparable with the values of other samples computed in this paper. The NEW 
and Burkert profile fits yield po = 2.04 x 10® M 0 kpc“® and 1.67 x 10® M 0 kpc“® 
respectively [16], thus confirming ( 21 ). 

To approximately quantify ppPP*^*^ in a familiar situation, we consider our 
own Milky Way galaxy as an example. Note that different groups have come up 
with somewhat different conclusions regarding the local (solar neighbourhood) 
density of dark matter [8-10]. The density, consistent with other standard esti¬ 
mates, seems to be Pq®'*^'' = (0.3 ± 0.1) GeV/cm® = 7.92 x 10® M 0 kpc“® [10]. 

that of K are not actually required as the behaviour of V" is practically insensitive to their 
variations within the Newtonian approximation. 
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Bergstrom, Ullio and Buckley [12] find local dark matter densities acceptable in 
the range 0.2 — 0.8 GeV/cm^, so that the fitted values roughly bunch around 
pSoiar ^ 2 ^q 6 M 0 kpc“^. For updated reviews, see [12,13]. With regard to the 
size of dark matter around our galaxy, note that the observed velocity profile is 
declining with radius: ’’The radial velocity dispersion shows an almost constant 
value of 120 km/s out to 30 kpc and then continuously declines down to 50 
km/s at about 120 kpc” [15]. Assuming a continuous fall, we can approximately 
take the dark matter radius to be Rdm ~ 180 kpc. Eq.(20) then yields a value 
PgPPer = 2.37 X 10^^ M 0 kpc“^. The actually fitted latest data on central density 
are = 4.13 x 10^ M 0 kpc“^ (Burkert profile) and 1.40 x 10^ M 0 kpc“^ 

(NFW profile) [16]. In both cases, we see that ppPP®*^ is about four orders of 
magnitude larger than the above fitted values of pa, confirming the inequality 
( 21 ). 

We should emphasize that it is the Doppler shifted light coming only from 
material stable circular orbits that led to the observation of anomalous rotation 
velocities [19]. That is precisely the reason why we focused on the circular orbits 
obtaining a stability inspired constraint. It is possible however that non-circular 
orbits or different kinds of instabilities may lead to stronger constraints than 
the ones derived in this work. For example it could be that galaxies with stable 
circular orbits, such as the Milky Way, actually present other unstable orbits 
and thus becomes non-viable within the EiBI scenario^] All kinds of orbits are 
of course dictated by the potential and in the present case the exact form, from 
Eq.(12), is 


^ (Po sin r - r) (r - po sin r -b ppr cos r) _ 


_9 

Lor 


r — Pq sin r p^r cos r 


ujr^ 


+E^ 


-1-f 


~ Po sin r)(r — p^ sinr + p^r cosr) 




( 22 ) 


where q = (i?DM/?’opt)^j = (-Rdm/t’")^, Po galactic sample 

characteristics, while and Lf are constants describing arbitrary particle tra¬ 
jectories. Once a sample is chosen, one can specify ranges for and Lf to see 
if the potential allows turning points so that non-circular orbits can lie between 
them (see, e.g., [20]). From the above, it is seen that the potential V{r) is a 
combination of periodic (sinf) and aperiodic (fcosr) functions. It turns out 
that, depending on the sample as well as on the prescription of and L^, 
the spacetime may or may not support turning points. Also, Eq.(ll) is non- 
integrable for this complicated potential so that orbit equations in closed form 
cannot be obtained. We believe that definitive conclusions about the possibility 
of other (un) stable orbits can nonetheless be drawn by adopting the more pow¬ 
erful method of autonomous dynamical system for Eq.(ll) with this potential. 
That would by itself be a separate task, which we leave for the future. 

In conclusion, we can say that the stability of orbits is a powerful constraint 
governing the dynamics of particles in the galactic halo (see e.g., [2,17-19]). 

®We thank an anonymous reviewer for pointing out this possibility. 
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The obtained stability inspired constraint (21) is a new testable prediction of 
the EiBI theory of dark matter developed by Harko et al. [4], The theory is 
increasingly gaining importance after it was originally developed by Bahados 
and Ferreira [21] on the basis of modified matter-gravity coupling [22-26]. It is 
powerful enough to handle different physical scenarios such as neutron stars and 
other compact objects first studied in Refs. [27-32]. Consistency with well known 
density profiles, as shown above, indicates that the EiBI theory could provide yet 
another viable alternative model for dark matter. However, a solution beyond 
the current Newtonian approximation and conclusively observed data on i?DM 
are the two ingredients that will be needed for a more precise testing of the 
theory in future. Also, note that if for just one galaxy, or cluster, the central 
dark matter density is higher than the predicted value ~ 10 ^^ M 0 kpc“^, 

then the EiBI scenario would fail as a candidate for explaining the anomalous 
rotation velocity of galaxies. 
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Figure captions 



Figure 1: Plot of V" vs R specific to the sample UGC 0128. The chosen 
parameters are: tq = 0.044, v^o = 0.000001, Topt = 27.6 kpc, k = 1.57 x 10^ 
kpc^ and pg = O.SOtt. The orbits are unstable in the chosen entire radial range 
R G [0.57r,7r] because V" > 0 there 
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Figure 2: Plot of V" vs R specific to the sample UGC 0128 The chosen param¬ 
eters are: tq = 0.044, Voo = 0.000001, Topt = 27.6 kpc and k = 1.57 x 10^ kpc^. 
Here central density is further lowered to Pq = 0.347r. The orbit is unstable in 
some intermediate radii as V" is partly positive and partly negative 



Figure 3: Plot of V” vs R specific to the sample UGC 0128. The chosen 
parameters are: tq = 0.044, Voo = 0.000001, ropt = 27.6 kpc and k = 1.57 x 10^ 
kpc^. Here central density is further lowered to Pq = 0.2997r, which corresponds 
to = 0.94 = /3“PP®’'7r, where /3 '^pp®’' = 0.299. The orbit is stable in the 

entire chosen range for R G [0.57r,7r]. 
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